Minimal model for beta relaxation in viscous liquids 
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Contrasts between beta relaxation in equilibrium viscous liquids and glasses are rationalized in 
terms of a double-well potential model with structure-dependent asymmetry, assuming structure 
is described by a single order parameter. The model is tested for tripropylene glycol where it 
accounts for the hysteresis of the dielectric beta loss peak frequency and magnitude during cooling 
and reheating through the glass transition. 

PACS numbers: 64.70.Pf, 77.22. Gm 



Viscous liquids approaching the calorimetric glass 
transition have extremely long relaxation times Q. The 
main relaxation is termed the alpha relaxation. There 
is usually an additional minor "beta" process at higher 
frequencies. Dielectric relaxation is a standard method 
for probing liquid dynamics ,2| . The study of dielectric 
beta relaxation in simple viscous liquids was pioneered 
by Johari and Goldstein more than 30 years ago 
but the origin of beta relaxation is still disputed 0,| 
It is unknown whether every molecule contributes to the 
relaxation |8J or only those within "islands of mobility" 
[fll llftllll. S imilarly, it is not known whether small angle 
jumps [8l. ll2l.ll3j or large angle jumps 0] are responsible 
for the beta process. 

Improvements of experimental techniques have re- 
cently lead to several new findings. The suggestion 
UH IHJ E3| that the excess wing of the alpha relax- 
ation usually found at high frequencies is due to an un- 
derlying low-frequency beta process was confirmed by 
long time annealing experiments by Lunkenheimer and 
coworkers [l8j ] (an alternative view is that the wing is 
a non-beta type relaxation process This lead to 

a simple picture of the alpha process: Once the effect 
of interfering beta relaxation is eliminated, alpha relax- 
ation obeys time-temperature superposition with a high 
frequency loss oc w^ 1 / 2 20j. Moreover, it now appears 
likely that all liquids have one or more beta relaxations 
[181 l2lL I22L |23j : Liquids like propylene carbonate, glyc- 
erol, salol, and toluene are now known to possess beta re- 
laxation, while o-terphenyl, previously thought to have a 
beta process only in the glassy state, has one in the equi- 
librium liquid phase as well. Finally, it has been shown 
that beta relaxation in the equilibrium liquid does not be- 
have as expected by extrapolation from the glassy phase: 
In some cases the beta loss peak frequency is tempera- 
ture independent in the liquid phase (e.g., sorbitol 15]), 
in other cases it is very weakly temperature dependent. 
On the other hand, the beta relaxation strength always 
increases strongly with temperature in the liquid phase 
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peak frequency and maximum loss for tripropylene gly- 
col cooled through the glass transition and subsequently 
reheated j2||. In the glassy phase (at low temperatures) 
the loss peak frequency is strongly temperature depen- 
dent while the maximum loss varies little. On the other 
hand, the loss is strongly temperature dependent in the 
equilibrium liquid phase. Here we even see the loss peak 
frequency decreasing upon heating. How is one to un- 
derstand these findings? A clue is provided by Fig. 1(b) 
which shows the maximum loss as a function of time af- 
ter an "instantaneous" temperature step, i.e., instanta- 
neous on the time scale of structural (alpha) relaxation. 
This experiment utilizes a special-purpose setup with a 
cell consisting of two aluminum discs separated by three 
capton spacers (layer distance ~ 20/x, empty capacitance 
30 pF). One disc, where temperature is measured via an 
NTC resistor, is placed on a Peltier element. Less than 
six seconds after a 0.672 K temperature jump is initi- 
ated, temperature is stable within 1 mK. In this setup 
we measure at 10 kHz which is the loss peak frequency 
(changes of loss peak frequency lead only to second order 
corrections of e^,.). The sampling time is 2 s. Figure 
1(b) shows a very fast change of the maximum loss, fol- 
lowed by relaxation toward the equilibrium value taking 
place on the structural (alpha) relaxation time scale. The 
existence of an instantaneous increase of the loss clearly 
indicates a pronounced asymmetry of the relaxing entity 
Inspired by this fact we adopt the standard asymmet- 
ric double-well potential model (Fig. 2) with transitions 
between the two free energy minima described by rate 
theory @, HHH |H 111 In terms 01 the sma11 barrier 
U and the asymmetry A, loss peak frequency / max and 
maximum loss e max are given [2°J by 
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The contrasts between beta relaxation in liquid and 
glass are clear from Fig. 1(a) which shows beta loss 
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The prefactor /o is assumed to be structure and temper- 
ature independent while €q(T) oc l/T [2^] is assumed to 
be structure independent: €q(T) = Tq/T. The free en- 



2 



4.06 - □ ft 



max 
^ £max 



(a) 

start, 7 end 



- -0.87 



N 

X 



3.96 - 



- -0.97 



J 3.86- 
J 3.76 - 



3.66 - 



end 



start 



- -1.07 &. 



-1.17 



--1.27 



bO 

o 

-J 



174 178 182 186 190 194 

T[K] 

i 

■1.18! 



■1.194 



r a 

^ -1.200 

so 
o 



■1.206 



■1.212 



AT=0.672 K 




— t<0;T=184.0K 

i i i i i 

1.00 3.00 5.00 

Log 10 (t [s]) 



FIG. 1: Data for tripropylene glycol, (a) Observation of beta 
loss peak frequency (□) and loss peak maximum (0) for a 
continuous passage through the glass transition starting from 
the equilibrium viscous liquid at 192 K, cooling to 174 K, and 
reheating to 194 K. After each temperature step the systems 
is kept at constant temperature for 50 min whereafter a spec- 
trum is measured (sampling time: 6 min) . In the glassy phase 
one observes the well-known strongly temperature-dependent 
loss peak frequency and weakly temperature-dependent max- 
imum loss; in the liquid phase the situation is the opposite 
(as seen also in other liquids, e.g., sorbitol [TH l24j|1. There 
is even a reversal so that the loss peak frequency decreases 
with increasing temperature, (b) Beta loss peak maximum 
monitored after a temperature increase of 0.672 K starting 
from equilibrium at 184.0 K. After 6 s temperature is stable 
within 1 mK of the final temperature. 



ergy differences U and A are expected to change with 
changing structure, but freeze at the glass transition. In 
terms of the fictive temperature Tf our model is based 
on 
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FIG. 2: Asymmetric double well potential as a simple model 
for beta relaxation. The two free energy differences U and 
A vary as structure changes with temperature in the equilib- 
rium liquid phase, but freeze at the glass transition. Working 
from this picture the simplest possible assumptions are that 
1) Structure is parameterized by just one order parameter, 
and 2) first order Taylor expansions apply in the relatively 
narrow temperature range studied. 



u 

A 



Uq + ak B Tf 
A - bk B T f 



(2) 



Equation (J2J follows from minimal assumptions: Suppose 
structure is parameterized by just one variable, s. Only a 
rather narrow range of temperatures is involved in studies 
of beta relaxation in the liquid phase and around the 
glass transition. Consequently, structure varies only little 
and, e.g., U(s) may be expanded to first order: U(s) = 
Cq+Ci(s — sq). For the equilibrium liquid, s — s(T) which 
may also be expanded to first order. By redefining s via 
a linear transformation we obtain s = T at equilibrium 
while U(s) is still linear in s. A single variable describing 
structure, which at equilibrium is equal to t emp erature, 
is - consistent with Tool's 1946 definition |3l| - to be 
identified with the fictive temperature: s = Tf. Thus 
one is lead to Eq. (J2J) where signs are chosen simply to 
ensure a, b > in fit to data. 

The model has 6 parameters: /o, To, Uq, Ao, a, and b. 
These were determined by measuring the instantaneous 
changes of f max and e^ ax upon a temperature step, as 
well as their thermal equilibrium changes. Figure 1(b) 
allows determination of the instantaneous change of the 
loss. It is not possible to determine the instantaneous 
change of / max - Instead we extrapolated measurements 
obtained by the standard cell |2fJ as follows (Fig. 3): 
Beta loss is monitored by first annealing at 183.0 K, sub- 
sequently changing temperature to 181.0 K. The latter 
data show a linear relation between log f max and log ej^ ax 
which, knowing the instantaneous change per Kelvin of 
e max from Fig. 1(b), is extrapolated to short times. 
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FIG. 3: Reversible temperature jump experiment for 
tripropylene glycol monitoring beta loss peak frequency and 
loss peak maximum. Starting at 185.0 K (not shown) temper- 
ature is first lowered to 183.0 K and kept there for 60 h; the 
last point, a, was obtained after further 24 h. Then temper- 
ature is changed to 181.0 K where it is kept constant for 140 
h (initially 20 min between measurements, later 12 hours). 
The point b is found by extrapolating to zero time after the 
quench utilizing the data of Fig. 1(b). This series ends at 
point c, thereafter temperature is changed back to 183.0 K 
where it is kept for 140 h (all points except the final point e 
refer to times up to 70 h; e was measured after further 70 h). 
Equilibrium at 183.0 K is somewhere between points a and e. 



In the data analysis it is convenient to eliminate cosh 
by introducing the variable 
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On the fast time scale U and A are frozen so Eq. (JTJ 
implies 
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Having determined the instantaneous changes of e: 
and / m ax, Eqs. and 0} provide 4 equations for the 6 
model parameters. The last two equations come from the 
temperature dependence of loss and loss peak frequency 
at thermal equilibrium where U (T) and A(T) are given 
by Eq. © with T f = T, leading to 
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Using Eq. 0} for the data of Figs. 1(b) and 3, and 
using Eqs. Q and Q for equilibrium state measure- 
ments at 183 and 185 K, the 6 parameters determined for 
tripropylene glycol H2 are: f = 6.2 • 10 11 Hz, T = 130 
K, U /k B = -1107 K, A /k B = 3039 K, a = 23.28, 
b = 12.77. Uq < reflects the fact that the beta loss 
peak frequency in the liquid phase decreases as temper- 
ature increases [Fig. 1(a)]. Physically, this anomalous 
behavior is caused by the barrier U increasing more than 
T upon heating. 

Once all parameters are fixed the model predicts 
how /max and e^ ax correlate for the continuous passage 
through the glass transition of Fig. 1(a). To analyze 
these data within the model we first note that if y(x) = 
l/cosh 2 (a;) then x = (f>(y) where <f)(y) = \a{l/^/y + 
\J\/y - 1), so Eq. (JTJ may be inverted: A/2k B T = 
0(e max /ef)'(T)) = Since lnY = -(2U + A)/k B T both 
<I> and lnY involve fictive temperature, and consequently 
both exhibit hysteresis at the glass transition. However, 
fictive temperature is eliminated by considering the fol- 
lowing variable 
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The model is tested in Fig. 4. The glass transition is not 
visible and most hysteresis is eliminated (better elimina- 
tion was obtained in Ref. [24| but without theoretical 
basis and with one free parameter [13|)- The line shown 
is the prediction of Eq. JH}. 

In conclusion, a minimal model for beta relaxation in 
viscous liquids has been proposed. The model is built 
on the four simplest possible assumptions: 1) Beta relax- 
ation involves only two levels, 2) structure is determined 
by just one order-parameter, 3) first order Taylor ex- 
pansions apply, 4) the two characteristic free energies U 
and A freeze at the glass transition. The model is clearly 
oversimplified. For instance, it predicts a Debye response 
which is not observed, and U and A would be expected 
to vary slightly with temperature in the glass. Never- 
theless, the model is able to rationalize the contrasts be- 
tween beta relaxation in liquids and in glasses. One final 
puzzling observation should be mentioned: The asym- 
metry A extrapolates to zero at a temperature which 
is close to the temperature where alpha and beta relax- 
ations merge. We have seen the same phenomenon in sor- 
bitol, a pyridine-toluene solution, polypropylene-glycol- 
425, and in 4,7,10-trioxydecane-l,13-diamine [lol |24| . 
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FIG. 4: Data for tripropylene glycol from Fig. 1(a) replotted 
to test the model. The variable Z defined in Eq. © is calcu- 
lated using parameters obtained by independent experiments. 



and have found no exceptions. This finding indicates that 
the merging temperature is fundamental, a symmetry is 
somehow broken below this temperature. 
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